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Abstract. The linear stability theory of Taylor-Couette flows (unbounded in z) is described includ- 
ing magnetic fields. Hall effect or a density stratification in order to prepare laboratory experiments 
to probe the stability of differential rotation in astrophysics. If an axial field is present then the mag- 
netorotational instability (MRI) is investigated also for small magnetic Prandtl numbers. For rotat- 
ing outer cylinder beyond the Rayleigh line characteristic minima are found for magnetic Reynolds 
number of the order of 10 and for Lundquist numbers of order 1. The inclusion of extra toroidal 
current-free fields leads to new oscillating solutions with rather low Reynolds numbers and Hart- 
mann numbers. The Hall effect establishes an unexpected 'shear-Hall instability' (SHI) where shear 
and magnetic field have opposite signs. In this case even rotation laws increasing outwards may be- 
come unstable. Recently global solutions have been found for the Taylor-Couette flows with density 
stratified in axial directions ('stratorotational instability', SRI). They exist beyond the Rayleigh line 
for Froude numbers of moderate order but also not for too flat rotation laws with fj^ < fl < fj. 



INTRODUCTION 

The rotation law in the infinite (in axial direction) Taylor-Couette flow is Q.{R) = 
a+b/R^, where a and b are related to the angular velocities ^2in and ^2out of the inner and 
the outer cylinders. The rotation laws with a = are called as located at the Rayleigh 
line. With i?in and 7?out as the radii of the cylinders the parameters fl = Qout/^^in and 
V = Rin/Roui of the flow are defined. According to the Rayleigh criterion the ideal flow 
is hydrodynamically stable when the specific angular momentum increases outwards, i.e. 
fi > f]^ or a > 0. The viscosity, however, stabilizes the flow so that it becomes unstable 
only if the inner cylinder rotates sufficiently fast. 

If the fluid is electrically conducting and an axial magnetic field is applied then after 
the old results the critical Reynolds number (for the inner rotating cylinder) grows with 
growing magnetic field. Figure 1 shows the theoretical results of Chandrasekhar (1961) 
for narrow gaps and very small magnetic Prandtl numbers, Pm = v/r] , together with the 
experimental data for Mercury of Donnelly & Ozima (1960). Theory and observations 
are in nearly perfect agreement with no indication of any magnetic-induced instability. 
The magnetic Prandtl number under laboratory conditions is really very small (10^^ for 
liquid sodium). 

Within the small-gap approximation but with free Pm Kurzweg (1963) found that for 
weak magnetic fields and sufficiently large magnetic Prandtl number the critical Taylor 
number becomes smaller than in the hydrodynamic case (Fig. 2). If the field is not too 
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FIGURE 1. Early results for small gap and small Pm (Chandrasekhar 1961). Left: The influence of the 
magnetic field (Q here the Hartmann number) on the critical Taylor number. Magnetic fields suppress the 
instabiUty. Right: The same for the wave number. The magnetic field elongates the Taylor vortices in the 
vertical direction. Boundary conditions: (a) conducting walls, (b) insulating waUs. 



Strong it can basically play a destabilizing role. 

Velikhov (1959) originally discovered this magnetic shear-flow instability which is 
now called 'magnetorotational instability' (MRI). He found that for the ideal hydromag- 
netic Taylor-Couette flow the Rayleigh criterion for stability changes to /i > 1 , i.e. only 
flows with superrotation are still stable. He found a growth rate along the Rayleigh line 
(i.e. a = 0) of IQi^f) and a critical wave number of 

k<2ftp, (1) 

with Va the Alfven velocity of the given axial field. Only if Va is smaller than the shear 
—R^dQ./dR the instability works which, therefore, is a weak-field instability. 

The hydrodynamic Taylor-Couette flow is stable if its angular momentum increases 
with radius, but the hydromagnetic Taylor-Couette flow is only stable if the angular ve- 
locity itself increases with radius. This remains true also for nonideal fluids. The MRI 
reduces the critical Reynolds number for weak magnetic field strengths for hydrodynam- 
ically unstable flow and it destabilizes the otherwise hydrodynamically stable flow for 
f)^ < ft < I. The MRI exists, however, in hydrodynamically unstable situations (jU < r)^) 
only if Pm is not very small (Fig. 2). 

As we shall demonstrate the magnetic Reynolds number 

Rm = Re Pm (2) 

controls the instability. Because of the high value of t] for liquid metals (exceeding 

1000 cm^/s), it is not easy to reach magnetic Reynolds numbers of the required order of 
1 ... 10. This is the basic reason why the MRI has not yet been observed experimentally 
in the laboratory. 
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FIGURE 2. Critical Taylor number as a function of the magnetic Prandtl number for strong and weak 
magnetic fields within the small-gap approximation. For weak fields and large magnetic Prandtl numbers 
the system is subcritical. For Bq — Q (or, what is the same, for Pm — 0) is Tacrk — 1750. (Kurzweg (1963). 




FIGURE 3. Bifurcation diagram for axisymmetric modes with resting outer cylinder of conducting 
material (left) and vacuum (right). Shaded areas denote subcritical excitation by the axial magnetic field. 
Tj = 0.5. From Rudiger, Schultz & Shalybkov (2003). 



MAGNETOROTATIONAL INSTABILITY (MRI) 

The Reynolds number is usually defined as Re = i?in(^out — ^in)^in/v. The ampli- 
tude of the external magnetic field Bq is expressed with the Hartmann number Ha = 
Bo\/ RiniRout — Rin) / (MoP^t]). The values of the Reynolds numbers above which the 
flow becomes unstable depend on the vertical wave number. They have a minimum at 
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FIGURE 4. Marginal stability lines for axisymmetric modes in containers with rotating outer cylinder 
of conducting material for Pm = 1 ( left) and Pm — 10^^ (right), fj = 0.5, /i =0.33. The instability domain 
is colored by grey and fluids in the cross-hatched area are always stable. 



some wave number for fixed other parameters. This minimum value is called the criti- 
cal Reynolds number. Figure 3 shows the neutral stability for axisymmetric modes for 
containers with both conducting and insulating walls with resting outer cylinder and for 
fluids of various magnetic Prandtl number. Re = 68 is the classical hydrodynamic solu- 
tion for resting outer cylinder and r) =0.5. There is a strong difference of the bifurcation 
lines for Pm > 1 (high conductivity) and Pm < 1 (low conductivity). For fluids with low 
electrical conductivity the magnetic field only suppresses the instability so that all the 
critical Reynolds numbers strongly exceed the value 68. 

For small magnetic Prandtl number the stability lines hardly differ. The opposite is 
true for Pm > 1. In Fig. 3 the resulting critical Reynolds numbers Re are smaller than 
68. The magnetic fields with small Hartmann numbers support the instability rather than 
to suppress it. This effect becomes more effective for increasing Pm but it vanishes for 
stronger magnetic fields. The MRI only exists for weak magnetic fields and high enough 
electrical conductivity and/or microscopic viscosity. 

Now the outer cylinder may rotate so fast that the rotation law no longer fulfills the 
Rayleigh criterion and a solution for Ha = does not exist. The nonmagnetic eigenvalue 
along the vertical axis moves to infinity but a minimum remains. Figure 4 presents 
the results for Pm = 1 and Pm = 10^. There are always minima of Re for certain 
Hartmann numbers (Riidiger & Shalybkov 2002). The minima and the critical Hartmann 
numbers increase for decreasing magnetic Prandtl numbers. For f\ = 0.5 and fl = 0.33 
the critical Reynolds numbers together with the critical Hartmann numbers are plotted in 
Fig. 5. Table 1 gives the exact coordinates of the absolute minima for experiments with 
rotating outer cylinder and for Pm = 10^^. They are characterized by magnetic Reynolds 
numbers of order 10, very similar to the values of the existing dynamo experiments. 



TABLE 1. Coordinates of the absolute minima in Fig. 4 for 
rotating outer cylinder with fi = 0.33, f) = 0.5 and Pm ~ 10^^ 
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FIGURE 5. The critical Reynolds numbers vs. magnetic Prandtl numbers marked with those Hartmann 
numbers where the Reynolds number is minimum, fj = 0.5, fi = 0.33. 



A container with an outer radius of 22 cm (and an inner radius of 1 1 cm) filled with 
liquid sodium requires a rotation of about 19 Hz in order to find the MRI. The required 
magnetic field is about 1400 Gauss. 

The results for containers with conducting walls are also given in Table 1. The 
minimal Reynolds numbers are higher than for insulating cylinder walls. The influence 
of the boundary conditions is thus not small^. 

There is a particular scaling for the special case of a = 0, i.e. for jl = i)^. One finds 
that the quantities ur^u^^Br and scale as VmT^/^ while a^,^^,^ and Ha scale as Pm^. 
The Reynolds number for the axisymmetric modes scales as Re Pm~^/^ (Willis & 
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Barenghi 2002). The scaling does not depend on the boundary conditions. However, for 
a > the much steeper scaling Re oc Pm^ results , leading to the surprisingly simple 
relation 

Rm ~ const. (3) 

for the magnetic Reynolds number Rm (Fig. 5). For the Lundquist number S = Hay/Pm 
of the characteristic minima we find 

S ~ const. (4) 

For small magnetic Prandtl numbers the exact value of the microscopic viscosity is thus 
not relevant for the instability. In consequence, the corresponding Reynolds numbers for 
the MRI seem to differ by 2 orders of magnitude, i.e. 10^ and 10^ at and shortly beyond 
the Rayleigh line. Figure 6 shows the behavior close to /i = f)^. There is a vertical jump 
within an extremely small interval of /i. This sharp transition does not exist for Pm = 1; 
it only exists for very small values of I*m. Even the smallest deviation from the condition 
jU = f)^ would drastically change the excitation condition. 
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FIGURE 6. Critical Reynolds numbers for the Taylor-Couette flow vs. p. for fj = 0.27 and Pm = 1 (left) 
and Pm = 10^^ (right). From all f), fj = 0.27 yields the lowest minimum. The curve for the hydrodynamic 
instability (Ha = 0) is dashed and the hydromagnetic curve (Ha > 0) is soUd. The dotted lines denote the 
location of a = 0. 

For f} ~ 0.5 the Velikhov condition (1) can also be written as kRi^ ~ Rm/S. The 
vertical extension 5z of a Taylor vortex is given by 



(5) 



The dimensionless vertical wave numbers k associated with the critical Reynolds num- 
bers are given in Fig. 7. For hydrodynamically unstable flows we have dz ^ Roat — Rm 
for small magnetic fields (Ha ~ 0). The cell therefore has the same vertical extent as it 
has in radius. 
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FIGURE 7. Wave numbers for which the Reynolds number is minimum for conducting walls and 
insulating walls, fj ~ 0.5, fl = 0.33, Pm = 10^^. Note how exact the Velikhov condition (1) is fulfilled. 
From Rudiger & Shalybkov (2002). 



The magnetic fields deforms the Taylor vortices. The deformation consists in a elon- 
gation of the cell in the vertical direction. The wave number is thus expected to become 
smaller and smaller for increasing magnetic field. This is indeed true for Pm 1, but for 
smaller Pm the situation is more complicated. 

The cell size is minimum for the critical Reynolds number for all calculated examples 
of hydrodynamically stable flows with a conducting boundary. This is not true, however, 
for containers with insulating walls, for which the cell size grows with increasing 
magnetic field. For experiments with the critical Reynolds numbers the vertical cell 
size is generally 2-3 times larger than the radial one. The smaller the magnetic Prandtl 
number the longer are the cells in the vertical direction. The influence of boundary 
conditions on the cell size disappears for sufficiently wide gaps. 

The given solutions of marginal stability of axisymmetric modes are stationary. One 
could ask for the character of the solutions if an extra toroidal magnetic field is applied 
(HoUerbach & Rudiger 2005). We know that current-free toroidal fields alone do not 
change the stability of Taylor-Couette flow (Velikhov 1959). If axial fields and current- 
free toroidal fields of the same order exist in the container, however, completely new 
solutions appear and they are oscillating (Fig. 8). The real parts of the eigenfrequencies 
are positive so that stationary modes do not longer exist. The Reynolds number of the 
flow and the Hartmann number of the axial field which are necessary to induce the MRI 
instability are strongly reduced by the toroidal field. In Table 2 characteristic numbers 
from numerical experiments with small magnetic Prandtl number and for fl = 0.3 are 
given. 




FIGURE 8. Critical Reynolds numbers for marginal stability of axisymmetric modes in containers 
(77 =0.5) with purely axial magnetic field (upper curves) and with extra (current-free) toroidal field 
(lower curves) for Pm = 10^^ and Pm = 10^^ as a function of p.. Note the extremely weak Prandtl 
number dependence of the lower curves and the strong reduction of the critical Reynolds numbers. From 
Hollerbach & Rudiger (2005). 



TABLE 2. Characteristic values for 
Tj = 0.5, A = 0.3, Pm = 10-5. J3 
describes the amplitude ratio of the 
toroidal and the axial magnetic field. 



Re Ha k 5K(co) 



1.6-10^ 850 1.7 
5.8-105 383 0.47 0.04 



SHEAR-HALL INSTABILITY (SHI) 

In the presence of the Hall effect there is even the possibility that differential rotation 
with Q. increasing outwards (i.e. dQ/dR > 0) becomes unstable. The Hall effect desta- 
bilizes flows with jl > 1 for which so far no other instability is known. The dispersion 
relation derived from the equation system with an induction equation including Hall ef- 
fect leads to the instability condition Rb dQ./ dR <0 with Rb as the ratio of the Ohmic 
time scale and the Hall time scale. It linearly depends on the electrical conductivity and 
the magnetic field (see Rudiger & Hollerbach 2004). For magnetic fields parallel to the 
rotation axis it might be positive and v. v. The result depends on the sign of the Hall 
resistivity; here the positive Hall resistivity is used. 

Figure 9 illustrates the instability for a container with again f) =0.5. The flow with 
jU < 1 is unstable for positive Rb. Flows with jU > 1 are unstable for negative Rb, 
i.e. if angular velocity and magnetic field have opposite directions^. The fact that the 
Hall effect destabilizes flows with the angular velocity increasing outwards was first 



For negative Hall resistivity the orientation is opposite 




FIGURE 9. The shear-Hall instability. Only the induction equation is solved for Taylor-Couette flows 
under the presence of the Hall effect. Left: magnetic field positive, shear negative. Right: magnetic field 
negative, shear positive. The Unes of marginal instabiUty are labeled with the actual value of A- 



described by Balbus and Terquem (2001). 

Figure 9 may only serve as an illustration. In particular for jl < 1 the situation is more 
complicated as simultaneously the MRI is active. One finds more details in the original 
paper of Rudiger & Shalybkov (2004) about the overall impact of the Hall effect on the 
Taylor-Couette flow. That the sign of the magnetic field has such an important influence 
on stability or instability of flows (e.g. Kepler flows) must have important consequences 
for astrophysical applications. 




FIGURE 10. The lines of marginal instabihty for p. =2 with axisymmetric (m = 0, solid line) and 
nonaxisymmetric (m = 1, dotted Une) modes. Left: conducting walls, note the crossover of both lines. 
Right: insulating walls, no crossovers. Here Ha = Rb. 

It is also important to know the nonaxisymmetric modes. After the Cowling theorem 
only nonaxisymmetric modes can be maintained by a dynamo process. We have dis- 
cussed the appearance of nonaxisymmetric modes for the magnetic Taylor-Couette flow 
with negative shear (Shalybkov, Rudiger & Schultz 2002). The common result was that 



the lines of marginal stability for m = and m = 1 have a very different behavior for 
different electrical boundary conditions. One finds crossovers of the stability lines for 
containers with conducting cylinder walls, and one never finds crossovers for containers 
with vacuum boundary conditions. The same happens for the Hall instability for mag- 
netic Taylor-Couette flows with /i > 1. In Fig. 10 the lines for both axisymmetric and 
nonaxisymmetric modes are given for conducting boundary conditions and for vacuum 
boundary conditions. Indeed, a crossover of the lines only exists for conducting cylinder 
walls. As usual, in the minimum the m = mode dominates but for stronger magnetic 
fields the mode with m = 1 dominates. 



STRATOROTATIONAL INSTABILITY (SRI) 

The hydrodynamic Taylor-Couette flow stability is now discussed under the presence of 
vertical density stratifications. Taylor-Couette flows with a stable axial density stratifi- 
cation were often studied with the conclusion that a stable stratification stabilizes the 
flow. Withjack & Chen (1974), however, reported new experiments with the Rayleigh 
line at /i = 0.04. For a slight negative density gradient disturbances of spiral form 
were observed for fl = 0.073, i.e. clearly beyond the Rayleigh line. Recently, Yavneh, 
Mc Williams & Molemaker (2001) find that the sufficient condition for stability against 
nonaxisymmetric disturbances is the same as the condition for MRI. Their numerical 
results demonstrated that certain instabilities exist also for finite viscosity. 




FIGURE 11. Taylor-Couette flow geometry with axial density stratification 

One has to look for the basic state with prescribed velocity profile and given density 
vertical stratification p = p{z). The initial system then can only exist if 



dz 



dR 



0. 



(6) 



is fulfilled. Under the centrifugal force the purely vertical stratification at the beginning 
transforms to a mixed vertical and radial stratification. This behavior strongly compli- 
cates the problem. 



In real experiments the initial vertical stratification is so small that R^Q. <^ g. Then 
after (6) the radial stratification is also small. Let us, therefore, consider the 'small- 
stratification case' with p = po + Pi(^,z) with pi <^ po, where po is the uniform 
background density and condition (6) is fulfilled at zero order. If only the terms of the 
zero order are considered then the system takes the Boussinesq form. is the vertical 
Brunt- Vaisala frequency A^^^ = —{g/ po)dpi/ dz and Fr = ^in/N is the Froude number. 

Here we suppose dpi /dz = const so that the coefficients of the system only depend on 
the radial coordinate and a normal mode expansion exp{i{m^ +kz+ cot)) is possible. As 
usual, only those wave numbers are considered which belong to the smallest Reynolds 
numbers. We have also checked the existence of the transition from stable to unstable 
state for some arbitrary points. 
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FIGURE 12. Left: The marginal stability lines for m > 0. The vertical lines mark the Rayleigh line 
fl — fj^ (left) and the outer limit line fi = i) (right). Right: The vertical wave numbers are very large. 
i) = 0.78, Fr = 0.5. From Shalybkov & Rudiger (2005). 

The dependence of the critical Reynolds numbers on ji is given in Fig. 12 (left). 
The already known marginal stability line for m = and also the new lines for m > are 
plotted. Note that the critical Reynolds numbers are remarkable small. The axisymmetric 
disturbances are unstable only for fi < i)^ in accordance to the Rayleigh condition. The 
new nonaxi symmetric disturbances are also unstable in the interval i)^ < fl < fj. The 
existence of the outer limit /i = r) is a new and surprising finding. Even more flat rotation 
laws with /i > f) are still stable (Shalybkov & Rudiger 2005). 

For m > the vertical wave number only slightly depends on m (Fig. 12, right). The 
wave numbers are so large that the cells are flat in the axial direction. The vertical wave 
number does hardly depend on fl. 

By Brandenburg & Rudiger (2005) the global problem is treated in Cartesian geome- 
try with the nonlinear Pencil code. For a Reynolds number of about 500 a resolution of 
128^ mesh points was sufficient. Again fj = 0.5. As expected, no instability is found for 
jU = 0.6. For smaller values { jl = 0.2) the flow is Rayleigh unstable. For fl = 0.3 the 
flow is indeed stratorotationally unstable for nonaxisymmetric modes also in the non- 
linear regime. The modes are saturated and they are azimuthally drifting. Cross sections 
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FIGURE 13. The global sti-atorotational instability for fi = 0.3, Fr = 0.32, Re = 500 and Pr = 10. Note 
the m = 2 azimuthal structure. From Brandenburg & Riidiger (2005). 



of Uz are shown in Fig. 13 for fl = 0.3 and Fr = 0.32. Note that there are 4 eddies in 
z-direction. The growth rate, however, with O.Olf^in proves to be rather small (see also 
DubruUe et al. 2005). 
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